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Abstract. We introduce a new model for the dynamics of a wormlike chain in 
an environment that gives rise to a rough free energy landscape, which we baptise 
the glassy wormlike chain. It is obtained from the common wormlike chain by an 
exponential stretching of the relaxation spectrum of its long-wavelength eigenmodes, 
controlled by a single parameter £ . Predictions for pertinent observables such 
as the dynamic structure factor and the microrhcological susceptibility exhibit the 
characteristics of soft glassy rheology and compare favourably with experimental data 
for reconstituted cytoskeletal networks and live cells. We speculate about the possible 
microscopic origin of the stretching, implications for the nonlinear rheology, and the 
potential physiological significance of our results. 
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1. Introduction 

Studies of in vitro polymerised networks and solutions of the biopolymers that constitute 
the cytoskeleton have provided many important insights into the molecular origin of 
the fascinating mechanical properties of cells and tissues [1]. Cell rheological data 
[2,3] are therefore often interpreted in terms of viscoelastic models developed for in 
vitro reconstituted cytoskeletal networks [4-8], and the inferred parameters have been 
suggested as pertinent indices for clinical diagnosis [9]. Also current rheological models 
for the active processes in the cytoskeleton are rooted in this viscoelastic paradigm 
[10-12]. In contrast, microrheological measurements [13] have revealed that cells obey 
the highly universal and comparatively featureless pattern of soft glassy rheology [14], 
ubiquitous in soft condensed matter [15]. Integrating these competing paradigms into a 
unified framework has become a major challenge in cell biophysics [16-18]. Recent 
experiments with highly purified reconstituted actin solutions have revealed strong 
signatures of a glass transition that might hold the key to a resolution of the dilemma 
[19]. The glass transition manifests itself in a strong stretching of the relaxation 
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spectrum, so that experiments probing the mechanical properties at a fixed time scale 
detect a sharp transition from fluid-like to solid-like behaviour as a function of various 
control parameters. In the following, we introduce a new model motivated by this 
observation, which is able to explain the apparently conflicting phenomenology on a 
common basis. We call it the glassy wormlike chain (GWLC), because its essence is 
an exponential stretching of the relaxation spectrum of the wormlike chain (WLC) 
model, which is the minimal model of a semiflexible polymer. The stretching is 
quantified by a single parameter £ , the stretching parameter, which can be thought 
of as a characteristic scale for the free energy barriers retarding the relaxation of the 
test chain's long wavelength eigenmodes. The characteristics of soft glassy rheology 
naturally ensue, including the widely reported small apparent power-law exponents 
in the frequency-dependence of the microrheological moduli [13,14,16,17], or the so- 
called "noise temperatures". Despite of the striking simplicity of the model {£ is the 
only free parameter), which is in stark contrast to prevailing viscoelastic models for 
cytosekeletal networks in vitro and in vivo [4,20], its predictions compare favourably 
with experimental data for live cells and reconstituted cytoskeletal networks. In 
particular, the logarithmic tails of the dynamic structure factor obtained by high- 
precision quasi-elastic light scattering from F-actin solutions at low temperatures and 
high concentrations [19] provide very direct evidence for the postulated exponential 
stretching. While a microscopic derivation of S remains so far elusive, we can 
demonstrate the practical usefulness of the model for many applications by computing 
pertinent measurable quantities. 

2. The wormlike chain (WLC) 

In the WLC model a semiflexible polymer is represented as a continuous space curve 
r(s) = (r_|_(s),s — T||(s)) with arc length s = . . . L. We consider the weakly-bending 
rod limit where deflections r' ± (s) from the straight ground state are considered small, 
r' ± (s) <C 1, the prime denoting an arc length derivative. From the arc length constraint, 
r /2 (s) = 1, it follows that longitudinal fluctuations are of higher order, r'» = 0(r' 2 ± ). 
The dynamics of a weakly bending WLC subject to an (optional) constant backbone 
tension / is to leading order described by the linear Langevin equation for its transverse 
excursions r±(s,t) 



Here k, (± and £±(s,t) denote the bending rigidity, the solvent friction per length, and 
Gaussian thermal noise, respectively. For more details about WLC dynamics the reader 
is referred to Ref. [21]. Equation (1) is solved by introducing eigenmodes, 



C±r L = -KrT" + /r'i + £ L . 



(1) 



oo 




(2) 



n=l 
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For simplicity we assume hinged boundary conditions, in which case the eigenfunctions 
are simple sine functions 

W n (s) = y|sin(/c n s) (3) 

and the eigenvalues k n = nn/L can be parametrised by natural numbers n. The 
eigenmodes relax independently and exponentially 

(a n (t)a m (0)) = S nm {al) exp(-*/r„) . (4) 

The equilibrium mode amplitudes 



follow from equipartition and the mode relaxation time is given by 

r n = r L /(n 4 + n 2 f/f L ) (6) 

with the relaxation time tl = ClL a / W 4 and the Euler force fi = ktt 2 /L 2 of the longest 
mode setting the characteristic time and force scale, respectively. 

These results can be used to calculate various time- dependent correlation functions 
as a superposition of eigenmode contributions. For example, the transverse dynamic 
mean-square displacement (MSD) reads 

SrUt) B (M, t) - r ,(s, „))») = g E 1 . (7, 

The MSD is directly or indirectly measured by a couple of experimental techniques, 
especially by particle tracking, dynamic light scattering, and diverse passive and (linear) 
active microrheology methods. 

3. The glassy wormlike chain (GWLC) 

The GWLC model is obtained from the WLC by an exponential stretching of the 
relaxation spectrum in the spirit of so-called hierarchically constrained dynamical 
models [22]. The strategy is also reminiscent of the generic trap models [23] underlying 
soft glassy rheology [14], but concerns the equilibrium dynamics of the test chain, here. 
The GWLC is a WLC with the relaxation times for all its eigenmodes of mode number 
n < I = L/A — or, more intuitively, of (half) wavelength \ n = L/n > A — modified 
according to 

~ / T n (n>l) 

n n I r n eMN n S) (n<l) ' W 

Here 

N n = l/n - 1 = A n /A - 1 (9) 

can be thought of as the number of interactions per length \ n with the environment and 
A <^ £ p , L as a characteristic interaction length. We moreover introduce the suggestive 
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notation ta = T\ = (^A 4 //t7r 4 and U\ = 2ti/t\ for the corresponding crossover time 
and frequency, respectively. We imagine the retardation of the relaxation of the long- 
wavelength modes of the test polymer to be caused by a complex environment such as 
a solution of other polymers. For example, in a semidilute solution of semiflexible 
polymers, A would correspond to the familiar entanglement length [24], u\ to the 
entanglement frequency, and N n to the number of entanglements an undulation of 
arc length \ n has to overcome in order to relax. With regard to the application 
to cytoskeletal networks we have in particular polymers with incompletely screened 
sticky interactions in mind (see below). The parameter £, which we call the stretching 
parameter, controls the slowing down caused by the interactions. It is the key parameter 
of the model — and in fact the only parameter apart from the rather obvious interaction 
scale A. Physically, it may be interpreted as a characteristic height of the free energy 
barriers in units of thermal energy k B T in a rough free energy landscape. Below, we 
speculate about the microscopic origin of such a free energy landscape in cytoskeletal 
networks and try to give tentative estimates for various contributions to £. 

Our definition of the GWLC does only affect the relaxation times but not the 
amplitudes a n of the test chain's eigenmodes. The question, how the amplitudes of the 
long wavelength modes are affected by the interactions with the disordered environment 
is in principle an interesting open question that is currently under investigation. 
However, we expect that the slowing-down captures the most important mechanism 
underlying the observed glass transition and the corresponding soft glassy rheology, and 
that the conformational aspects are in this respect of minor relevance. 

As a first example of a pertinent observable for the GWLC, we plot in figure 1 (left) 
the MSD for various values of the stretching parameter at vanishing prestress. 

4. Prestress 

Before we come to the evaluation of further observables, we first want to consider the 
effect of tension on a GWLC. An (optional) constant tension / was already included 
in our brief account of the ordinary WLC in section 2. However, we should certainly 
also expect an effect of any kind of external or internal stress onto the escape of the 
polymer over the free energy barriers represented by £. Intuitively, the force is expected 
to "help the polymer over the free energy barriers", but we cannot, of course, exclude 
the opposite effect, namely that the traps become under certain circumstances deeper 
upon applying a force. In any case, the natural way to introduce a force into this picture 
is via a tilting of the free energy landscape in the spirit of a Kramers escape rate model, 
i.e. 

5-£±/// T , f T = k B T/A. (10) 
The minus sign corresponds to the force lowering the barrier. The length A should be 
interpreted as a characteristic width of the free energy wells and barriers. Accordingly, 
fx represents the scale of thermally induced force fluctuations, which are even present 
in absence of an applied stress. 
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Figure 1. Time dependence of the transverse dynamic mean-square displacement 
8r\ L (t) (MSD) and the dynamic structure factor S(q,t). Left: The normalised 
MSD of a GWLC as a function of time t for stretching parameters £ — 5, 15, 35 
(from top to bottom) for / = 0. The straight line indicates the limiting power-law 
growth 5r\ L (t) oc t 3/4 of the MSD of an infinite weakly bending WLC. Right: The 
dynamic structure factor of a GWLC, evaluated numerically for q 2 = € p 7r 4 /A 3 and 
stretching parameters £ = 0,5,15,35 (from left to right). The straight lines indicate 
the logarithmic intermediate asymptotics calculated in equation (13). 

The introduction of an external force in equation (10) may be seen as a simple 
heuristic method to address the rheology of prestressed networks or even of the 
nonlinear rheology of cytoskeletal networks. This is of considerable interest for potential 
applications of our model, since prestress is thought to be the crucial element needed for 
mimicking typical cell rheological behaviour using much simpler reconstituted networks 
[18, 25-28]. There is, however, a subtle point involved in the interpretation of the force / 
as arising from a prestress. For simplicity, we have tacitly assumed above that the force 
that pulls the test chain over the free energy barriers is identical to its backbone tension. 
This should not a priori be a problematic assumption for qualitative purposes. But a 
test polymer embedded into a cytoskeletal network or sticky biopolymer solution that 
was initially unstressed will generally have a conformation different from the stretched 
conformation of our free test chain equilibrated under the backbone tension /. As noted 
above, we do not attempt to address this problem, at the present stage. Altogether, after 
disregarding such potential modifications of the equilibrium conformation of the test 
chain by its surroundings and identifying the barrier lowering force with the backbone 
tension, / remains to be related to the macroscopic shear stress. Consistent with our 
discussion of the shear modulus in section 7, below, we follow Ref. [29] in writing 
a = //5£ 2 , where £ = (3/c p L) 1//2 is the mesh size of a semidilute solution of semiflexible 
polymers. Its relation to the polymer concentration c p and the numerical prefactors are 
of geometric origin. As a reminder of the tentative nature of the identification of o as 
an actual prestress we call a the "nominal prestress". 
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5. Dynamic structure factor 



For sufficiently long times t ^> (±/ nq 4 , the dynamic structure factor S(q, t) at scattering 
vector q 3> A -1 follows from the transverse mean square displacement according to [30], 

S(q,t)/S(q,0) ~ exp [- q V ±L (t)/4] . (11) 

The time dependence of the structure factor is plotted in figure 1. A pronounced 
logarithmic intermediate asymptotics is seen to develop for large S. It can be calculated 
approximately as follows. First, we identify a logarithmic intermediate asymptotics in 
the MSD. At intermediate times the MSD takes the form 

4A 3 

±LW~^i A (oo) + ^^[ 7E + log(VTA)] (l«t/r A «£), (12) 



Sr 2 



£i p 7l 4 



for an unstressed chain (/ = 0). Here <5r^ A (oo) contains the saturated contributions 
from the free WLC modes up to wavelength A. The GWLC part of the mode spectrum 
is asymptotically approximated by the logarithmic term. Expanding the exponential in 
equation (11) to leading order in the logarithmic contribution gives 



S(q,t) 
S(q,0) 



2A3 



q 2 A 
££ p ir 4 



7s + log — 

TA 



exp 



q 2 5r 



_LA 



OO 



(13) 



for the unstressed test chain. As demonstrated in figure 1, the numerically evaluated 
dynamic structure factor agrees with this approximation well beyond the time domain 
where the logarithmic contribution in the exponent of equation (11) is small. From 
the slope of the logarithmic tails of the structure factor in a semi-logarithmic plot the 
stretching parameter S is thus immediately inferred. 



6. Microrheology 

From the transverse MSD of a point on the polymer contour, we deduce the linear 
susceptibility oif(u) (the subscript / refers to the prestressing tension) to a transverse 
oscillating point force at frequency uj from the fluctuation dissipation theorem. It relates 
the imaginary part &"f(uj) of the susceptibility to the Fourier transform 8r\ L {uj) of the 
MSD via ctf(u) = —uj5r 2 _ L (uj)/2kBT. The real part of aj is then uniquely determined 
by the Kramers-Kronig relations, so that we find altogether 

af(uj) = fc^V* £ (n 4 + n*f/f L )(l + i^fn) ' (W) 

For better comparison with the macrorheological complex shear modulus, it is customary 
to report the inverse (up to a constant scale factor) g*f{uj) oc l/aif(u) of the susceptibility, 
which is called the "mi crorheo logical modulus". Its real and imaginary parts g'f(ui) and 
g"f{u>) are plotted in figure 2. The prestressing force / is seen to compete with the 
stretching parameter S in raising/lowering the apparent power-law exponent of the low- 
frequency modulus. Its full effect is somewhat richer, because / also affects the WLC 
mode amplitudes and relaxation times according to the explicit expressions provided in 
section 2. 
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Figure 2. Frequency dependence of the real and imaginary parts g'Aui), g'j(ui) 
(red/blue curve) of the microrheological modulus, which is a common representation 
of the dynamic linear response to a transverse point force applied to a test polymer. 
The abscissa is normalised to 5j =0 (^a) (£ = 0). Left: £ = 3 (solid), 10 (dashed), 
35 (dotted) at vanishing prestress, / = 0. Right: f / /a = (solid), 1 (dashed), 2.4 
(dotted) at fixed £ — 25 with fx — 0.1/a, /a being the Euler force for the interaction 
wavelength A. 



7. The shear modulus 

The microrheological modulus discussed in the preceding section should not be confused 
with the shear modulus measured by a macroscopic rheometer. While, in practice, the 
two quantities are sometimes hard to distinguish, this must be attributed to non-ideal 
(i.e. not point-like) probes such as colloidal beads used to transmit the force to the 
medium and to detect its deformation, which require additional considerations [31,32]. 
Noninvasive techniques, such as quasi-elastic light scattering are more sensitive to the 
difference between the two response functions [19]. In the following, we present results 
based on the assumption that the macroscopic shear modulus is obtained by applying 
our GWLC prescription, equation (8), to the high frequency limiting form of the shear 
modulus [29] . The latter is a single polymer quantity due to the independent relaxation 
of the short wavelength modes that dominate the high frequency response. The results 
thus obtained for the frequency dependence of the real and imaginary parts G'^iu) 
and G"(u) of the shear modulus G*(u) are displayed for various nominal prestresses 
a = //5£ 2 in figure 3. In contrast to the microrheological modulus, the shear modulus 
is seen to develop a plateau that is sensitive to the prestress as a consequence of the 
underlying assumption [29] that the single polymers are stretched affinely upon applying 
a macroscopic shear stress. This assumption is presently under scrutiny [33-35], and 
might in the future have to be relaxed for modes of wavelength A n > A. 

On the level of our simplifying assumptions, the shear modulus in the presence 
of a prestress is equivalent to the nonlinear differential shear modulus K [5,36]. It is 
therefore of some interest to evaluate the dependence of \G*\ at a fixed frequency u as 
a function of the nominal prestress a, which then amounts to a characterisation of the 
nonlinear finite-time elasticity of the system. Via the force dependence of the bare WLC 
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Figure 3. Frequency dependence of the real and imaginary parts G^(w), G^(w) 
(red/blue curve) of the macroscopic shear modulus obtained by applying the GWLC 
prescription to the theoretical expression for the high-frequency limiting form of the 
shear modulus from Ref. [29] for various nominal prestresses a = //5£ 2 corresponding 
to /// A = (solid), 2 (dashed), 2.45 (dotted) at £ = 25, f T = 0.1/ A . The modulus 
has been normalised to G' a=0 (iOA) (£ = 0). 



expressions of section 2, the prestress causes stress stiffening, i.e., a monotonic increase 
of \G*\ with a. In contrast, the exponential speed-up of the relaxation caused by the 
barrier-lowering effect of the generated tension according to equation (10) eventually 
overcompensates this stiffening for large stresses. This is analysed in figure 4 (left) 
for various stretching parameters. The limiting functional relation \G* C \ oc cr 3//2 for the 
stiffening, which is only slowly approached for £ — > oo is an immediate consequence 
of the underlying affine assumption [29]. While there is some recent experimental 
support that this limiting behaviour is indeed measurable in actin bundle networks 
heavily crosslinked by scruin [36] (plausibly corresponding to 8 — > oo), experiments 
for actin/a— actinin solutions [37], actin solutions homogeneously crosslinked by heavy 
meromyosion (HMM) in the rigor state [8], and pure actin solutions [19] rather reveal a 
continuity of stiffening relations. They suggest that the stiffening is much less universal 
than previously thought [6] and might depend on the degree of "glassiness" of the 
system, consistent with a finite stretching parameter £ dependent on various control 
parameters such as crosslinker concentration, temperature and ionic strength. To some 
extent, observations of a weaker stiffening might also indicate a contribution of (non- 
affine) transverse modes, as these exhibit weaker stiffening; see the corresponding curves 
for the transverse microrheological response in figure 4 (right), which converge to the 
asymptotic stiffening relation \g*j\ oc / for very large £. But the GWLC does neither 
require nor support a strict correspondence of network affinity to the nonlinear response 
as recently postulated [5]. 



The Glassy Wormlike Chain 



9 




Figure 4. The normalised moduli as a function of prestress. Left: The shear modulus 
|G* |/|G*_ |, which is, under the assumptions stated in the main text, equivalent to 
the nonlinear differential shear modulus, evaluated at a fixed frequency to = O.Iwa 
as a function of the normalised tension ///a corresponding to a nominal prestress 
a = //5£ 2 for various £ = 4. . .40 (from bottom to top) at fa = 0.1 /a- The straight 
line indicates the asymptotic stiffening power-law to which the stiffening curves slowly 
converge in the limit £ — > oo. Right: The corresponding curves for transverse the 
microrhcological modulus. 

8. Cytoskeletal networks and rough free energy landscapes 

While free energy landscapes certainly represent an appealing intuitive framework for 
the discussion of many properties of complex disordered media, it is notoriously hard 
to calculate their pertinent features from first principles [38]. In the remainder, we 
want to give some qualitative arguments as to why we expect the rough free energy 
landscapes alluded to in our motivation of the GWLC to be characteristic of cytoskeletal 
networks in vitro and in vivo. We distinguish two contributions to the free energy: direct 
contributions from a bare polymer-polymer pair interaction potential, and collective 
cageing or entanglement effects. For the typical biopolymer solutions and networks 
found in cells and tissues both cannot easily be derived microscopically or reduced to 
known elementary atomic pair interactions: the former because of the strongly collective, 
heterogeneous, anisotropic etc. nature of protein interactions [39], the latter because 
semidilute polymer solutions represent a highly correlated state of matter [40]. 

While protein interactions remain poorly understood, many observations hint at 
(unspecific, e.g. hydrophobic) adhesive contact interactions incompletely screened 
by electrostatic repulsion [41]. It is thus not implausible that direct interactions of 
cytoskeletal elements can approximately be represented by a pair potential U(r) of 
the qualitative form sketched in figure 5 [42], which features a narrow energy barrier. 
Such "enthalpic" barriers slow down the mode relaxation by an Arrhenius factor, which 
scales exponentially in the barrier height. The barriers would thus yield substantial 
contributions to £ without seriously affecting the thermodynamics of the system. The 
less obvious free energy contributions from cageing and entanglement are typically 
postulated in generic free- volume theories [43], as their calculation from first principles 
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Figure 5. Schematic sketch of a hypothetical bare pair potential U(r) suggested to 
provide an approximate representation of the complex interactions between cytoskeletal 
constituents. 

remains one of the major challenges in the theory of structural glasses [44]. For 
semiflexible polymer networks they can be estimated within the tube model [45] , which 
suggests a contribution £ ~ 1. Accordingly, a plausible estimate for the effective well 
width A should be given by the tube diameter, which scales like £, 6 ^ 5 £p ^ 5 [24] in a 
semidilute semiflexible polymer solution and is estimated to assume values in the range 
of some 10 to 100 nm for typical in vitro polymerised actin solutions. In any case, 
simple exponential scaling of the relaxation times in £ and in the wavelength A, as 
postulated in equation (8) seems plausible for an entangled solution of weakly bending 
sticky polymers. 

It is important to realise that the parameters £ and A are to be understood as 
effective parameters. They cannot generally be expected to correspond directly to 
some salient features of a bare interaction potential as sketched in figure 5 nor to the 
experimentally more accessible coarse-grained interactions between adjacent polymer 
segments of length A. (The latter can be roughly thought of as a smeared-out version 
of the former, due to the effect of thermal undulations of the polymers and possible 
compliant crosslinkers [46].) 

The above suggestion to interpret £ essentially as a kind of kinetic "stickiness" 
parameter for cytoskeletal polymers might have interesting implications as to the 
physiological role played by the broad class of actin binding proteins such as crosslinkers 
and molecular motors. While little is known about the microscopic interactions between 
cytoskeletal polymers, even less is known about how they might be affected by these 
proteins. Yet it seems plausible that the predominant effect of most crosslinkers can be 
subsumed into the parameter £. Also, while current approaches to the active rheology of 
cells [10, 11,47] put much emphasis on the dynamic role of molecular motors in causing 
transport and motion, some of them might well play a less flamboyant role most of the 
time: namely adjusting the effective stickiness £ and tension o to tune the viscoelastic 
properties of the cytoskeleton while saving it from glassy arrest [48]. 
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9. Conclusion 

In summary, we have presented a simple modification of the standard model of a 
semiflexible polymer (the WLC), which we call the glassy wormlike chain (GWLC), 
and which is obtained by an exponential stretching of the relaxation spectrum of the 
WLC. The striking resemblance of the predicted mechanical observables with rheological 
data for cells and reconstituted cytoskeletal model systems naturally suggests that 
these systems must exhibit strong rheological redundancy. In fact, the perplexing 
universality and robustness of the mechanical performance of biological cells and tissues 
against structural modifications by drugs and mutations has been an enigma in cell 
biology for quite some time [13,49]. Our glassy wormlike chain model offers a very 
economical explanation in terms of the difference between the stretching parameter £ 
and the prestress a (in natural units). If changes in the structure and interactions of 
the cytoskeleton affect its rheological properties chiefly via £ and a, the potential for 
mutual compensation appears indeed enormous. This leads us to the suggestion that a 
microscopic calculation of the various anticipated contributions to £ will be one of the 
major challenges for future theoretical work aiming to explain the rheology of cells and 
in vitro reconstituted models of the cytoskeleton. 
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